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Abstract 

Coupled spin-charge drift-diffusion equations are derived for a biased two-dimensional electron 
gas with weak Rashba spin-orbit interaction. The basic equations formally agree with recent results 
obtained for spin-orbit coupled small polarons. It is shown that effects of an in-plane electric field 
on a homogeneous spin system can completely be described by an associated in-plane magnetic 
field. Exploiting this analogy we predict among other things the electric- field equivalent of the 
Hanle effect. 

PACS numbers: 72.25.-b, 72.10.-d, 72.15.Gd 
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I. INTRODUCTION 



The prospects of a new generation of electronic devices stimulates a renewed theoretical 
and experimental interest in the study of spin effects in semiconductors. A prerequisite for 
the design of semiconductor structures, whose function is based on electron spin, is a due 
understanding of spin dynamics and spin-polarized transport. In this field, the spin-orbit 
interaction (SOI) received particular interest since it allows purely electric manipulation of 
the electron spin. Many studies refer to a two-dimensional electron gas (2DEG), in which the 
Rashba SOI arises because of the quantum well asymmetry in the perpendicular direction. 
The SOI leads to a coupling between spin and charge degrees of freedom, which offers the 
possibility of controlling the spin polarization by an electric field. 

The theoretical description of spin phenomena in semiconductors under the influence 
of SOI is based on appropriate transport equations. There are numerous approaches 
that rely on a separation of drift-diffusion processes for spin and charge 
densities. Strictly speaking, such an approach is inappropriate, when the SOI has to be ac- 
counted for. Spin-charge coupled drift-diffusion transport equations have been derived from 
firm microscopic models for extended states of a 2DEG la, Is, 10, 11] and for the hopping 



transport of small polarons. 12| Especially, the approach by Mishchen 



so et al 



14 



m 



for a 



16] Unfortunately, 



2DEG with Rashba SOI initiated many interesting studies, 
variants of these calculations suffer from an inconsistency because SOI contributions from 
the collision integral have been disregarded. The corrections lead to a cancellation of con- 
tributions that has been used in a previous work [13j to erroneously predict propagating 
coupled spin-charge waves. 

The aim of the present paper is to derive spin-charge coupled drift-diffusion equations for 
a 2DEG with Rashba SOI that correct deficiencies of previous approaches and that describe 
the influence of an external in-plane electric field on spin polarization. With the help of 
universal macroscopic drift-diffusion equations, the spatial and temporal evolution of coupled 
spin-charge disturbances as well as associated charge accumulation and magnetization are 
studied in semiconductor heterostructures with Rashba SOI. Our basic equations, which are 
derived for weak SOI, completely agree with results obtained for the hopping transport of 
small polarons. [12] The field-induced homogeneous spin accumulation as well as the charge- 



Hall current are treated. Furthermore, it is shown that for a homogeneous system, the effect 
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of the electric field on spin polarization can be completely captured by a fictitious magnetic 
field. This analogy between the real applied electric field and an auxiliary magnetic field is 
used to predict a number of interesting electric-field effects on spin. We mention the decay 
of a spin polarization by a transverse electric field. Using this electric-field driven Hanle 
effect by exchanging the in-plane magnetic field by an electric field in the measurement set 
up for the ordinary Hanle effect, one should be able to alternatively determine electron and 
spin lifetimes under steady-state conditions by varying the electric field strength. Another 
application refers to the pseudo charge-Hall effect [17|, which is induced by circular polarized 
light via the creation of a permanent spin magnetization. 

II. KINETIC EQUATIONS 

The effect of an in-plane electric field E on coupled spin-charge excitations of semi- 
conducting electrons in an asymmetric quantum well can be described by a single-particle 
Hamiltonian 

H = 22 a lx [ £ h ~ £f] OfcA - ^ (fojfc • ctav) a[ x a ky 

k,X k,X,\' 

2 « | K= 

fc.A k,k> X 

which includes both the Rashba SOI and the short-range spin- independent elastic scatter- 
ing on impurities. The Hamiltonian is expressed by creation (a kX ) and annihilation (g^a) 
operators that depend on the vector k = (k x , k y , 0) and the spin index A. We introduced 
the Fermi energy ep, the vector of Pauli matrices er, and the strength u of the short-range 
elastic impurity scattering, which is characterized by the momentum relaxation time r. The 
energy e k of free electrons in the 2DEG and the coupling term of the Rashba SOI are given 
by 

h 2 k 2 J~i , , . ma . . 

Sk = -z — , uj k = —(Kxk), K=— e z , (2) 
2m m nr 

where m denotes the effective mass and a the strength of the SOI. The central quantity of 
our approach is the spin-density matrix /, the components of which 

ti(k,k , \t) = (ai x a k , y ) t , (3) 

satisfy kinetic equations, which are derived for the corresponding physical elements / = Tr/ 
and / = Trer/. The four components of the spin-density matrix depend on two wave vectors 



(k + k')/2 — > fc and (fc — fc') — > k, which allow the description of inhomogeneous charge and 
spin distributions. Relying on the Born approximation for the treatment of elastic impurity 
scattering and restricting to lowest-order corrections of the SOI to the collision integral, we 



obtain the following Laplace-transformed kinetic equations 



id 



sf - -(k • fc)/ - -K(f x «) + 6 -^V k f = -(f -f ) + fo, (4) 

mm n t 

sf + 2K x /) - -( K • fc)/ + —(K x K )f + ^-V k f 

mm n 

r rd£ fc r cte fc 

with the initial charge and spin distribution fo and fo, respectively. An integration over 
the polar angle (p of the vector fc is indicated by a cross line over respective quantities and 
s denotes the Laplace variable that refers to the time t. The second and third term on the 
right-hand side of Eq. (jSJ) stems from spin contributions of the collision integral, which have 
to be taken into account to guarantee that the spin system correctly approaches the state 
of thermodynamic equilibrium. The solution of the coupled integro-differential Eqs. (jl]) and 
flSJ) is searched for in the long-wavelength and low-frequency regime. 

To illustrate our approach, let us first treat the evolution of charge- density disturbances 
at zero SOI. As the inelastic scattering time t £ is usually much larger than r, the quasi- 
momenta thermalize already at the time scale r E > t > r (sr 1). During this period, the 
density matrix (/, /) approaches its mean value (/, /) with respect to the angle if. In the 
following stage of the evolution t > r e , which lasts till a characteristic diffusion time r^, the 
carrier density locally approaches the equilibrium distribution. Consequently, the behavior 
of particles in this time interval can be described by the Fermi function n(Ek) with a Fermi 
energy that depends on spatial coordinates r and time t: f(£k, r I t) = n(ek — ep(r,t)). In 
this evolution period, the energy is already thermalized although both the charge and spin 
densities remain still inhomogeneous. For weak perturbations Ep(r,t) = Ep + Aep(r, t), we 
obtain for the density fluctuation 5/(e fc ,r | t) = —AsF(r,t)dn(Ek — Ep)/dEk- This result 
brings us to the separation ansatz 

for the new unknown function F(n \ s), where n = J dep(e)n(e) with p{e) being the density 
of states of the 2DEG. For brevity, we write e instead of and use a prime to indicate a 



derivative with respect to e. It is in line with this discussion and Eq. ([6]) to replace the drift 
term accordingly 

^E x ^-f^eE x ^cos(<p)J. (7) 
h ok x m n 

Adopting these approximations, which express the basic understanding of the drift-diffusion 
approach, Eq. (BJ is easily solved under the condition of vanishing SOI (K = 0). A spectral 
drift-diffusion equation is obtained by expanding the solution of Eq. (j3J) with respect to n 
and by integrating over the angle (p 



n" 



s + D(k)K 2 + iue — E ■ k 



/(e,«|s)=/o, (8) 



TV 

where the diffusion coefficient and the mobility are given by D(k) = (Hk) 2 T/(2m 2 ) and 
H = er/m, respectively. The final integration over the energy e leads to the well-known 
drift-diffusion equation 

[ a + Dk 2 - ifjiE ■ k] /(« | s) = / , (9) 

for the charge density f(n \ s) = J dep(e)f(e, k \ s). The relationship between the diffusion 
coefficient D and the mobility [i is given by the Einstein relation fi = (eD jn)dnj 'dep, which 
is applicable both for Fermi and Boltzmann statistics. 

Within the framework of the drift-diffusion approach, a similar approximation can be 
used for the spin components of the density matrix, when the Rashba SOI is weak (Q = 
<^ 1, which is accessible by tuning the SOI coupling constant a via the shape of the 
confinement potential). In this case, the spin relaxation time r s is large so that we can focus 
on the time hierarchy r < r £ < t s < r^. Under the condition t > r e but t/r s arbitrary, a 
nonequilibrium spin polarization exists on the background of thermalized carrier energies. 
Therefore, Eqs. ([6]) and (|7j) can be used also for the spin contributions in the kinetic Eqs. (jl]) 
and (jSJ). Exploiting these approximations, the following set of linear equations is obtained 
for the components of the spin-density matrix 

r , r^f J- r \ 2eEr Tl" . 

af + iU(q x fy - q y f x ) = R ^T £ ~^7 cos (v)/ 

2 n" — 2eEr n" — 
af x + 2tt cos(y?)/ z - iQq y f = R x + -e— sm(<p)f — e— cos{ip)f x 



7 n' hk n' 

2 n" , ^ 2eEr n" 

-e—cos((p)f — — e— ; 

7 n' hk TV 

2eEr n" 



af y + 2ttsm(tp)f z + iVLq x f = R y e— cos{<p)f — cos(ip)f y 



af z - 2n(cos((f)f x + sm(ip)f y ) = R z f^T e ~^ cos (^)A ? ( 10 ) 
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with q x>y = K XtV /k, 7 = k/K and 

a = <Jo — i^VL(q x cosy? + q y sin <p), cr = st + 1. 



(11) 



It is assumed that the in-plane electric field E is oriented along the x axis. Restricting to 
lowest-order contributions in n and E, we obtain 



R* = f s +Tf t 



xQ 



ihKr {en')' — 
—K y f } R = f + Tf , R z = f t 



VTLOq TV 



rf, 



- ihKr {en')' - f hKr Aen")'- 
f y + rf x0 + r — —K x f eE — /. 



(12) 



(13) 



mcr n' 

It is straightforward but cumbersome to solve the equations for /, /, to expand the solution 
with respect to k, and to calculate the final integral over the angle ip. What we obtain 
by this procedure are spectral drift-diffusion equations for coupled spin-charge excitations. 
After integrating over the remaining energy e, we get our final result 



d_ 
dt 



- i(n ■ fiE) + Dk 2 



F — —n ■ [K x F] + ^ (k ■ [K x fiE]) (K-F)= 0, (14) 

m e 



d A 
— -i{K- fiE) + Dk 2 + — 

dt t s 

+ ^^( K .[Kxf,E])K- 



F + AD 

2h/x 



K x [in + —E 

2D 



x F 



K x (in + —E) 

v 2D ' 



0. 



(15) 

< er s i 2L> J J 

with A xx = Ayy = 1, A zz = 2, and the spin-scattering time calculated from l/r s = ADK 2 . 
In Eq. fll5p . the electric field is accounted for via the quasi-chemical potential (in — > in + 
jiE/(2D)). These coupled spin-charge drift-diffusion equations are valid for weak SOI (Q -C 
1 so that t/t s 1). What is interesting is that Eqs. fTl4l) and ffl5|) formally agree with 
results, which were recently derived by a different approach for the hopping transport of 



small polarons. 



12] Spin effects in the latter system exclusively occur in the weak SOI 



regime as the lattice constant is much smaller than typical values of K . Summarizing this 
observation, we point out that the drift-diffusion equations have an universal character for 
the Rashba model with weak SOI. 



III. RESULTS AND DISCUSSION 



The afore-mentioned analogy with the hopping transport provides us a recipe to transfer 
results recently obtained for small polarons 12| in a straightforward manner to spin effects 

6 



of extended electronic states. We shall not repeat all calculations already presented in Ref. 
[12I but restrict ourselves to some additional conclusions. 

To begin with, let us treat the field-induced homogeneous (k = 0) spin accumulation. 
From Eq. (IT5|) . we obtain for the steady-state field-mediated magnetic moment 



/ = hfi(K x E)ri. 



(16) 



This result expresses the well-known magnetoelectric effect that was predicted by Edelstein 
[19I many years ago. For its derivation, it was essential to account for spin contributions on 
the right-hand side of Eq. (jSJ), the origin of which is the collision integral. If we would neglect 
these corrections, a term of the form ihKn x f jm remained uncompensated in Eq. (TT5]) . which 



was used in a recent paper 



131 ] to predict coupled spin-charge waves. Introducing an electric 



n 



field via the quasi-chemical potential [18j], this term would lead to a magnetic moment 
/ ~ 1/K, in complete disagreement with well established results l^ . 

There are numerous other field-mediated spin effects that follow from Eqs. (|14p and (1151) . 
As an example, we treat spin waves that exist in a stripe of a 2DEG oriented parallel to 
the electric field (and the x-axis). Taking into account t/t s <C 1 and restricting ourselves 
to long wavelengths (Dk^t s <C 1), we obtain the dispersion relation 



^1,2 



3 

2n 



% ± ^/Ts/te - 1 



i± 



TE 



(17) 



in which a rate 1/te — i^E) 2 / (9D) appears, which is associated with the electric field. Field- 
induced damped oscillations arise at sufficiently high electric field strengths (t s /te > 1). For 
typical parameters \i ~ 10 6 cm 2 /Vs, t s « 40 ps, r = 0.5 ps, and n = 10 11 /cm 2 , we obtain the 
condition E > 3.25 V/cm, which gives rise to an appreciable current density in the 2DEG. 
A much stronger electric field (t s /te 1) drives spin-charge coupled oscillations with the 
constant field-dependent frequency uj =\ fj,E \ / (2^Dt s ). In contrast to the well-known 
space-charge waves of free and trapped electrons, this massive mode with the frequency 
uj = KfiE is independent of the propagation vector k. 

To get further information on spin-charge coupling due to an electric field, the drift- 
diffusion Eqs. (fl4"]l and f|T5|) are expressed in spatial coordinates. For the charge density, we 
obtain the continuity equation 

dF(r,t) 



dt 



+ dwj(r,t) = 0, 
7 



(18) 



with a particle current 



h 2hu 2 
j = (uE - DV r )F + — \K x F] —K [K x E\ F z , 

m e 

that includes both charge and spin components. Besides the charge Hall current 



(19) 



j H = -2h» 2 K [KxE]F z 



(20) 



which arises from a given out-of-plane spin polarization F z , there appears another spin- 
related term that is responsible for the spin-galvanic effect. As the spin is not conserved, 
the equation for the spin components F a of the density matrix cannot be written in the form 
of a continuity equation. Rather, we obtain from Eq. (fT5|) 



+ —F a + 2/x {[K x E] x F) a - ^- [K x E] a F + ^ = 0, (21) 



dt t s 

where the spin current J ia is given by 



(fiE - DV r ) F a + 5 J a 



(22) 



with the spin components 



5J n = —(K x F) 



2hKr 



m 



m 



K x fiE)F g , SJ Z 



2Dm 



K x fiE)F 



5J X 



hK 2 r 



m 



y ) ' 



5J„ 



HK 2 t 



m 



F z Sy 



(23) 



Studies of the spin current received a great deal of recent interest in the literature. [18|, [20 , 



21 



22|, 



23|, 



24 



251 ] This interesting discussion is confronted with a serious problem, namely 



that the spin is not conserved. As a consequence, different definitions of the spin current 
have been put forward in the literature. According to Eq. (|23l) . we obtain for the spin-Hall 
current J z = —2hK 2 fi 2 EF/e, which is neither universal 26| nor in line with the result 



derived from a more physically motivated definition of the spin-Hall current. [18|, [20j In 
addition, there is a component of the spin current that is completely independent of the 
electric field (J® = —(hK/m)(T/T s )F). This astonishing result is replaced in an alternative 
approach 18| by a spin current that is related to the initial variation of the spin accumulation 
and that disappears, when the latter one reaches its steady-state value. As the concept of 
the spin current is not well founded, it seems to be expedient to avoid its introduction. 



S 



To proceed, we further focus on a homogeneous electron gas, for which the effect of an 
electric field on the spin polarization can completely be simulated by an appropriate in-plane 
magnetic field of the form 

H cS = ^[KxE). (24) 
Replacing the electric field by this equivalent magnetic field, Eq. (|2~T!) is written as 

OF A _ 2«b rrr ^ 1 xH e g „ 

— + -F + [H eS x F - _^L_!^ = o, 25 

with % = /i^n' being the magnetic susceptibility. The term xH e s is responsible for the 
spin accumulation, whereas the vector product H c g x F leads to spin precession around the 
effective magnetic field H c g. This close relationship between spin polarization due to an 
electric field and its description by an associated in-plane magnetic field can be used for the 
derivation and interpretation of electric-field effects on spin. As an example, we mention the 
rotation of an initial perpendicular homogeneous spin polarization F z (t = 0) = F z q into the 
plane of the 2DEG due to Larmor precession. For a constant electric field oriented along 
the x axis, we obtain from Eq. (I25I) 



with H c q = hfiKE j hb- This solution demonstrates that sufficiently high electric fields lead 
to an in-plane spin polarization that oscillates with the frequency uj s w 2fiBH e s/h = 2fiEK. 

To further exploit the analogy between electric and magnetic field effects, let us treat the 
optical generation and recombination of a steady-state spin polarization under the additional 
influence of a real in-plane magnetic field H oriented along the y axis. In this case, Eq. (|25|) 
takes the form 

^ + ^ F + ^£[ Bxi F] = G -^, SF = F-**, (27) 
at t s h r n B 

where B = H + H c g. The vector G describes the optical out-of-plane spin generation and 
To is the relaxation time of photogenerated electrons. The steady-state solution of Eq. (|27|) 
is easily obtained 

t±G z 11112 1 

with u c = 2fiBB/h. At zero electric field, the solution describes the depolarization of spin 
(and hence the degree of circular polarization of the luminescence) by a transverse magnetic 



field, which is known as the Hanle effect. 



271 ] Combining the measurement of the zero-field 



spin generation G z with the magnetic-field dependence in Eq. ( 128]) (Han 



e effect), both the 



electron lifetime and the spin- relaxation time t s can be determined. [28|, |29|, |30j, [31|, [32|, 133 ] 
According to Eq. §ZEi), an effect of the same kind exists also at zero magnetic field (H = 0) 
due to an in-plane electric field. To describe this effect, it is only necessary to replace 
the Larmor frequency 2fiBB/h by uo c = 2fiEK in Eq. (|28|) . This Hanle effect driven by 
a pure in-plane electric field can likewise be used to measure lifetimes of charge and spin 
excitations. To the best of our knowledge, this promising experimental technique has not 
been employed hitherto. With respect to the charge transport, there is the interesting effect 
that a circular polarized light illumination induces a pseudo-Hall effect 17j in the absence of 
any external magnetic field. A quantitative description of this Hall contribution is provided 
by Eqs. fl20|) and fl28l) . Another application of the electric- magnetic field correspondence 
refers to a modification of recent experiments, in which the optically induced spin-galvanic 



effect was measured. 



34j Instead of using an external magnetic field to achieve an in-plane 



spin polarization necessary for the occurrence of the effect, one can likewise apply an in-plane 
electric field. 



IV. SUMMARY 

Based on the density-matrix approach, spin-charge coupled drift-diffusion equations were 
derived for extended electronic states in a 2DEG with weak SOI. The final basic equations 
agree with results that were recently obtained for the hopping transport of spin-polarized 



polarons. [12J Due to this correspondence, results on spin transport obtained for localized 
and extended states are mutually applicable to each other. In the course of the derivation, it 
was found that a consistent treatment of spin-charge coupling requires a careful consideration 
of spin-orbit contributions to the collision integral, which give rise to a tricky cancellation 
in transport equations. Disregarding these corrections is the source of fatal mistakes that 
plague former approaches. 

Particular emphasis was put on the effect of an electric field on the spin polarization of 
a homogeneous 2DEG. It was shown that the electric field can be replaced by a fictitious 
in-plane magnetic field in the drift-diffusion equations for the spin components. This inter- 
pretation of the equations reveals a number of interesting similarities between spin effects 



10 



induced by electric or magnetic fields. From an experimental point of view, most attractive 
seems to be the electric-field equivalent of the Hanle effect, which provides a new possibility 
to measure lifetimes of spin and charge excitations by manipulating exclusively an in-plane 
external electric field. 
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